In-Medium Similarity Renormalization Group for Nuclei 
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We present a new ab initio method that uses similarity renormahzation group (SRG) techniques to 
continuously diagonalize nuclear many-body Hamiltonians. In contrast with applications of the SRG 
to two- and three-nucleon interactions in free space, we perform the SRG evolution "in medium" 
directly in the A-body system of interest. The in- medium approach has the advantage that one 
can approximately evolve 3, A-body operators using only two-body machinery based on normal- 
ordering techniques. The method is nonperturbative and can be tailored to problems ranging from 
the diagonalization of closed-shell nuclei to the construction of effective valence-shell Hamiltonians 
and operators. We present first results for the ground-state energies of ^He, ^^O and ^"^Ca, which 
have accuracies comparable to coupled-cluster calculations. 

PACS numbers: 21.30.Fe, 21.60.De, 21.10.-k, 21.60.Cs 



Great progress has been made in ab-initio nuclear 
structure over the past decade, where it is now possi- 
ble to calculate properties of light nuclei up to about 
carbon [1 , ^2,] and low- lying states of medium- mass nuclei 
near closed shells [sf. A key challenge in nuclear physics 
is to extend this ab-initio frontier to larger and open- 
shell systems. This requires methods that can handle 
the strong coupling between low and high momenta in 
nuclear forces used in these calculations. 

In recent years, new approaches to nuclear forces based 
on renormalization group (KG) ideas have been devel- 
oped that decouple high-momentum degrees of freedom 
by lowering the resolution (or a cutoff) scale in nuclear 
forces to typical nuclear structure momentum scales 
Such RG-evolved potentials, known generically as "low- 
momentum interactions," greatly simplify the nuclear 
many-body problem and enhance the convergence of 
structure and reaction calculations, while the freedom 
to vary the resolution scale provides a powerful tool to 
assess theoretical uncertainties due to truncations in the 
Hamiltonian and from many-body approximations 0-0 • 

One path to decouple high-momentum degrees of free- 
dom is the similarity renormalization group (SRG), which 
was introduced independently by Glazek and Wilson [H 
and Wegner [9|. The SRG consists of a continuous 
sequence of unitary transformations that suppress off- 
diagonal matrix elements, driving the Hamiltonian to- 
wards a band- or block-diagonal form. Writing the uni- 
tarily transformed Hamiltonian as 

H{s) = U{s)HU\s) = H'^is) + i/°^(5) , (1) 

where H^{s) and H'^^{s) are the appropriately defined 
"diagonal" and "off- diagonal" parts of the Hamiltonian, 
the evolution with the flow parameter s is given by 

^ = [^is),H{s)]. (2) 



Here r]{s) = [dU{s)/ds] U^s) is the anti-Hermitian gen- 
erator of the transformation. The choice of the generator 
first suggested by Wegner, 

vis) = [H''{s),H{s)] = [J?d(s), , (3) 

guarantees that the off-diagonal coupling of is driven 
exponentially to zero with increasing s [9|. Through dif- 
ferent choices for and one can tailor the SRG 
evolution to transform the initial Hamiltonian to a form 
that is most convenient for a particular problem 0, |Tl| . 
It is this flexibility, together with the fact that one never 
explicitly constructs and applies the unitary transforma- 
tion U (s) [rather it is implemented implicitly through the 
integration of Eq. (|2])] that makes the SRG a powerful 
alternative to conventional effective interaction methods 
such as Lee-Suzuki similarity transformations [4]. 

To date, the SRG applications to nuclear forces have 
been carried out in free space to construct "soft" nucleon- 
nucleon (NN) and three-nucleon (3N) interactions to be 
used as input in ab-initio calculations [3, ^] • While the 
free-space evolution is convenient, as it does not have to 
be performed for each different nucleus or nuclear mat- 
ter density, it is necessary to handle 3N (and possibly 
higher-body) interactions to be able to lower the cutoff 
significantly and maintain approximate cutoff indepen- 
dence of A ^ 3 observables. The SRG evolution of 3N 
operators represents a significant technical challenge that 
has only recently been solved in a convenient basis [6|. 

An interesting alternative is to perform the SRG evo- 
lution directly in the A-body system of interest (9l-[ll|. 
Unlike the free-space evolution, the in-medium SRG (IM- 
SRG) has the appealing feature that one can approx- 
imately evolve 3, ...,A-body operators using only two- 
body machinery. The key to this simplification is the use 
of normal-ordering with respect to a finite-density ref- 
erence state. Starting from a general second-quantized 
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FIG. 1: Convergence of the in- medium SRG results at the 
normal-ordered two-body level, IM-SRG (2), for ^He using 
the generators 77^ (left) and 77^^ (right panel). The filled 
(open) symbols correspond to solving Eqs. (|9|)- (pT]) with 
the underlined terms omitted (included). The ground-state 
energy ^o(oo) is given as a function of the harmonic os- 
cillator parameter hu with increasing single-particle space 
Gmax = max(2n + /). The initial NN interaction is a free-space 
SRG-evolved potential with A = 2.0 fm~^ from the N^LO po- 
tential of Ref. [13]. For comparison we show the coupled- 
cluster CCSD and CCSD(T) energies in the Cmax = 8 space 
(calculated at their huo minimum). 

Hamiltonian with two- and three-body interactions, all 
operators can be normal ordered with respect to a finite- 
density Fermi vacuum |<l>) (e.g., the Hartree-Fock ground 
state), as opposed to the zero-particle vacuum. Wick's 
theorem can then be used to exactly write H as 

H = Eo^Yl ^^i^^'} ^^Yl ^^^^^ {ala]aiak} 



ijkl 



X! ^iJkimn {ala^jala„amai} , 



3!2 



(4) 



ijklmn 



where the normal-ordered strings of creation and anni- 
hilation operators obey {^\{al ■ ■ ■ aj}]^) = 0, and the 
normal-ordered 0-, 1-, 2-, and 3-body terms are given by 



Eo 



fij 
^ijkl 



i ij 

I "^ifkijk , (5) 

ijk 

Tij+Y.^!:Sk^k + \Y.vil,inuni, (6) 

k kl 
y(2) , 1 T.(3) 



ijklmn 



(3) 



ijklmn ' 



(8) 



Here, the initial n-body interactions are denoted by V'^^^ 
and rii = 6{eF — £i) are occupation numbers in the ref- 



erence state 1^), with Fermi energy e^. It is evident 
from Eqs. (j5j)-([7]) that the normal-ordered terms, Eq^ f 
and F, include contributions from the three-body interac- 
tion V^^^ through sums over the occupied single-particle 
states in the reference state |<l>). Therefore, truncating 
the in-medium SRG equations to normal-ordered two- 
body operators, which we denote by IM-SRG (2), will ap- 
proximately evolve induced three- and higher-body inter- 
actions through the nucleus-dependent 0-, 1-, and 2-body 
terms. As a preview, we refer to Fig. [T] with the very 
promising convergence of the ^He ground- state energy, 
which is comparable to coupled-cluster results. 

Using Wick's theorem to evaluate Eq. (|2]) with H{s) = 

truncated to 



(1) 



V 



Eo{s) + f{s) + r(s) and rj = rj 
normal-ordered two-body operators, one obtains the cou- 
pled IM-SRG (2) flow equations (with = 1 — n^): 



dEo (1) ^ / \ I -'- (2) - - 

= l^Vij'fji {rii - Uj) + -2_^%ki^kiij nnfijUknu 

ij ijkl 

(9) 

E[^llV. + (i-2)] 

i 

f Y^^^ ~ ^^(^I}^ r^-ii2 - fij r]f^i2) 

ij 



+ ^ E [v'kiij^ijk2{ninjnk + niUjUk) + (1 o 2) , (10) 

ijk 



dT 



= E - fli r/Ll) -(1^2) 

y:[ fa^3^ri2.4 -fer?i4)-(3^4) 

i 

ij 

- "i) [(€i4rai3 - n%^j2iA) -(1^2) 



(11) 



The IM-SRG (2) equations exhibit important similarities 
to the CCSD approximation of coupled-cluster theory. 
For instance, the commutator form of the flow equations 
gives a fully connected structure in which H{s) has at 
least one contraction with rj. Therefore, there are no un- 
linked diagrams and the flow equations are size extensive. 
Combined with the 0{N^) scaling with the number of 
single-particle orbit als, this makes the method well suited 
for calculations of medium-mass nuclei. The IM-SRG 
is intrinsically nonperturbative, where the flow equa- 
tions, Eqs. (|9])- (pT]) , build up nonperturbative physics 
via the interference between the particle-particle and 
the two particle-hole channels for F and between the 
two-particle-one-hole and two-hole-one-particle channels 
for /. The perturbative analysis reveals that the IM- 
SRG (2) energy is third-order exact (as is the CCSD 
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FIG. 2: In- medium SRG evolution of normal-ordered two- 
body matrix elements Tijki connecting hole- hole (hh) and 
particle-particle (pp) states for ^^O starting from a smooth- 
cutoff Mow k with A = 1.8 fm~^. The color scale is in MeV, and 
initial and s-evolved results are shown. The axes label two- 
body jj-coupled states \{na,laja,tza), (ribjbjb.tz^); J = 0). 
The Tijki where ijkl = ppph or hhhp, which are not driven 
to zero with the current generator, are not shown. 

approximation) and that / and F are second-order ex- 
act It also implies that for calculations with harder 
interactions, the underlined terms in Eqs. ([9])- (fTT]) should 
be excluded because they produce higher-order contribu- 
tions (with alternating signs) to Eq that are also gen- 
erated by the inclusion of higher-body normal-ordered 
interactions, rj'^^^ and corresponding to simultane- 
ous 3p3h excitations. Because such triples excitations 
can be sizable for hard potentials, the underlined terms 
in Eqs. (|9])- (pT]) should be omitted to better preserve 
the partial cancellations that would occur against the 
[7^(3) contributions. This is consistent with the ob- 
servation in Fig. [1] that for soft potentials our results are 
insensitive to the inclusion of these terms. Therefore we 
define the IM-SRG(2) truncation without these terms for 
consistency. 

In this initial study, we restrict our attention to 
the ground states of doubly-magic nuclei and define 

^od(^) = /°^(5) + r°^(5), with 

ph 

T°'^{s)= rpp,hh'{s){alal,ah'ah} + li.c., (13) 

where p^p' and h^h' denote unoccupied (particle) and 
occupied (hole) Hartree-Fock orbit als, respectively. We 
consider two different cases for the generator rj. First, we 
take the Wegner choice r]^{s) = [i^^(s), i^°^(5)]. Second, 
we follow White jlllj and define 



ph 



fph {ala-h} 



fp-fh- r 



H.c. 



phph 



, I'pp^hh^ialalah'ah} 

pPhh' ^ - Jh- Jh' + ^pp'hh' 
where App/fiiii = Tppfpp/ + ^hh'hh' ~ ^phph ~ ^p'h'p'h' ~ 
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FIG. 3: Convergence of the IM-SRG(2) energy Eo{oo) for "^He 
using the generator 77^^ and starting from the "bare" N^LO 
potential of Ref. 13]. The notation is the same as in Fig. [T] 
The converged IM-SRG(2) energy agrees well with the CCSD 
result (the coupled-cluster energies are taken from Ref. [Sj]), 
while second- and third-order many-body perturbation the- 
ory, MBPT(2) and MBPT(3), clearly break down. 



^ph'ph' - ^p'hp'h and fp = fpp (the s dependence is sup- 
pressed for simplicity). Both generators suppress off- 
diagonal {Iplh and 2p2h) couplings and drive the Hamil- 
tonian towards diagonal form, 



H{oo) = ^o(oo) + /^(oo) + F^(oo) , 



(15) 



but White's generator (77^^) is significantly more efficient, 
because the flow equations are less stiff in this case and 
the evaluation of r] at each step is significantly faster. 
The evolved Hamiltonians using rf and rf^ are unitar- 
ily equivalent if no truncations are made. Any differ- 
ences in energy eigenvalues therefore provide a measure of 
the truncation error resulting from neglected three- and 
higher-body normal-ordered terms in our calculations. 

At the end of the flow, the reference state becomes 
the ground state of i^(oo), with fully interacting ground- 
state energy Eo{oo)^ and |$) decouples from the rest of 
the Hilbert space (Ipl/i, 2p2h^ . . . , ApAh sectors), 



QH{oo)P = and PH{oo)Q = 0, 



(16) 



where P = \^){^\ and Q = 1 - P. This decou- 
pling follows from the observation that all other normal- 
ordered couplings annihilate the reference state, [f^{s) -h 
F^(s)]|<l>) = 0. Combined with /°^(oo) andF°^(oo) being 
driven to zero, this implies the block-diagonal structure 
of Eq. (p!6|) . The IM-SRG is very flexible and alterna- 
tive choices of H'^^ (and r]) can be used to target excited 
states, single-particle properties, and to construct effec- 
tive valence shell-model Hamiltonians and operators for 
open-shell systems [lo|, Ell • 

Figure [1] shows the IM-SRG (2) ground-state energy 
Eo{oc) for ^He calculated in increasing spaces defined by 
the single-particle emax = max(2n + /). For all cases the 
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FIG. 4: Convergence of the IM-SRG(2) energy ^o(oo) for ^^O 
(left) and ^°Ca (right panel) using the generator 77^^ (solid 
symbols) and in comparison to CCSD results (dashed lines). 
The notation is the same as in Fig. [H The initial Vnn is 
a smooth-cutoff V^iowfe with A = 1.8 fm~^ for ^^O and a free- 
space SRG potential with A = 1.8 fm~^ for ^°Ca, both evolved 
from the N^LO potential of Ref. [H. 

flow equations, Eqs. (|9])- (pT]) , were solved in a jj-coupled 
basis. The r]^ and 77^^ results agree within 20 keV, which 
suggests the truncation to normal-ordered two-body in- 
teractions is a controlled approximation. This is consis- 
tent with Fermi system arguments for interparticle in- 
teractions where a finite-density reference state is close 
to the interacting ground state [15]. In addition, the 
IM-SRG(2) emax = 8 energy is essentially converged and 
within 20 keV of the exact NCSM diagonalization [6], and 
in good agreement with the coupled-cluster CCSD(T) en- 
ergies (based on the code of Ref. [3]). We stress that the 
agreement is obtained at the normal-ordered two-body 
level without including residual three-body interactions. 

The suppression of H^^{s) is illustrated in Fig. [21 
which shows the 77^-evolution of normal-ordered two-body 
matrix elements Tijj^i. As expected, the off-diagonal cou- 
plings {ijkl = pphh or hhpp) are rapidly driven to zero. 
An important practical consequence is that many-body 
approximations become more effective under the SRG 
evolution before complete decoupling has been reached. 

Figure [3] shows the IM-SRG(2) results for ^He starting 
from a "bare" N^LO potential, which is a harder initial 
interaction. The ground-state energy clearly converges to 
a value close to the CCSD result. The failure of many- 
body perturbation theory in this case verifies that the 
IM-SRG is an intrinsically nonperturbative method. 

Finally, we apply the IM-SRG to calculate the ground- 
state energies of ^^O and ^^Ca in Fig. HI As for the 
^He results of Fig. [3l the calculations are well converged 
and have accuracies that closely track the CCSD ener- 
gies. As discussed above, the IM-SRG (2) includes some 
simultaneous 3p3h excitations for Eo{s) that partially 
cancel against contributions that would arise if normal- 



ordered three-body operators were kept in the flow equa- 
tions. This motivated excluding the underlined terms in 
Eqs. ([9|)- (pT]) . The omitted terms are negligible for soft 
interactions, as shown in Fig. [H but they become larger 
for hard interactions such as the "bare" N^LO potential 
used here, and thus require a consistent treatment either 
by omitting them in the IM-SRG (2) equations, or by in- 
cluding normal-ordered three-body operators in the flow 
equations. In the former case, we find here an accuracy 
that is comparable to CCSD calculations. 

In summary, we have shown that the in-medium SRG is 
a promising method for ab-initio calculations of light and 
medium-mass nuclei. The use of normal ordering allowed 
us to evolve the dominant induced 3, ...,74-body interac- 
tions using only two-body machinery. We have presented 
first IM-SRG (2) results for the ground-state energies of 
closed-shell nuclei, which were in very good agreement 
with CC calculations. Work is in progress to include 3N 
forces and to study effective valence shell-model Hamil- 
tonians and operators for open-shell systems. 
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